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Abstract 

Relativistic exchange current corrections to the impulse approximation in low 
and intermediate energy neutrino-nucleus scattering are presented assuming 
non-vanishing strange quark form factors for constituent nucleons. Two-body 
exchange current operators which treat all SU (3) vector and axial currents 
on an equal footing are constructed by generalizing the soft-pion dominance 
method of Chemtob and Rho. For charged current reactions, exchange cur- 
rent corrections can reduce the impulse approximation results by 5 to 10 % 
depending on the nuclear density. A finite strange quark form factor may 
change the total cross section for neutral current scattering by 20% while ex- 
change current corrections are found to be sensitive to the nuclear density. 
Implications on the current LSND experiment to extract the strange quark 
axial form factor of the nucleon are discussed. 
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I. INTRODUCTION 



The existence of two-body meson exchange currents (MEC) in electromagnetic and weak 
axial interactions in nuclei has by now been firmly established [0-^ . In the electromagnetic 
sector the need for two-body currents has been realized even before the discovery of the pion 
by Siegert in 1937 through his consideration of vector current conservation [^. However, 
it took almost 35 years before the first solid quantitative evidence for MEC was discovered 
when the 10% discrepancy between theory and experiment in the threshold radiative np 
capture rates was remedied in terms of one-pion exchange current correction The ex- 
planation soon after of the cross section for the inverse reaction, the electrodisintegration 
of the deuteron, in terms of MEC [^] left little doubt of the important role that two-body 
currents play in electromagnetic interactions of the deuteron. For larger nuclei, it is well 
known that effects of MEC are best found in low and intermediate energy magnetic isovector 
processes. These two-body currents play important roles in realistic descriptions of diverse 
nuclear phenomena such as the renormalization of orbital (^-factors, 6gi [^, magnetic form 
factors of light p-shell nuclei [|], transverse (e, e') response functions in the dip region [^-[11 



and cross-sections for electromagnetically induced two-nucleon emission reactions |T^ . 

The role of MEC in weak axial transitions in nuclei has been investigated by Kubodera, 
Delorme and Rho who predicted a large renormalization of axial charge due to two-body 
currents. An ideal place to test their prediction is in the first -order forbidden /3-decays whose 
transition amplitude involves a cancellation between time and space components of the one- 
body nucleon current, thus making it sensitive to two-body MEC effects. Indeed, Guichon, 
GifFon and Samour [l^ found that the impulse approximation prediction of the ratio of 
/i-capture to /9-decay rates for the 0"^ 0~ transition between ^^O and ^^A^ was a factor of 
2 larger than the measured value. This large discrepancy disappeared when they included 
two-body MEC corrections to the impulse approximation. In addition, recent shell model 
analysis of first forbidden /3-decay rates covering a wide range of nuclei [[15HT7[| indicates 
a substantial exchange current contribution to the renormalization of weak axial charge in 
nuclear medium, thus confirming the prediction of Kubodera, Delorme and Rho. These 
solid empirical evidences of two-body currents in electromagnetic and weak axial processes 
in low and intermediate energy nuclear phenomena strongly motivate to investigate whether 
charged and neutral currents, where both vector and axial currents are involved, are also 
subject to renormalizations in nuclei due to MEC. 

Another reason to study MEC effects in low and intermediate energy neutrino-nucleus 
scattering is that it has been receiving increasing attention as a means to measure the 
strangeness matrix elements of the nucleon [lT8|-p5[] . The measurement of polarized structure 
function gi and the extraction of the sum rule suggested the possibility of a rather large 



strange quark axial matrix element for the proton leading to the so-called "spin crisis" p6 



Although there are numerous works attempting to understand the role of hidden flavor in 



nucleons [27|, the situation regarding the strangeness degrees of freedom in the nucleon is 



far from clear and it is hoped that neutral current neutrino-nucleus interactions might be 



able to shed a new light into this problem p8 



For example, Garvey et al. pointed out that the ratio of proton-to-neutron yield in 
quasi-elastic neutral current neutrino-nucleus scattering, hereafter denoted by R{p/n) = 
a{v, z/'p)/(t(z/, v'n), is an observable that is sensitive to the strange axial form factor of the 
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nucleoli p2| , |2^ . This ratio is currently being measured in the LSND experiment at Los 
Alamos [2^. In their work Garvey et al. calculated R{p/n) within a non-relativistic RPA 
framework using the impulse approximation and later included the effects of final state in- 



teractions experienced by the ejected nucleon with a continuum RPA formalism pj]. This 
correction was found to have about 30% effect on individual neutrino-nucleus cross sec- 
tions but cancelled out in the ratio of proton-to-neutron yields. A similar calculation by 
Horowitz et al. using a Relativistic Fermi Gas (RFG) model in the impulse approxima- 
tion did not include the final state interactions but the resulting R{p/n) was found to be 
similar to that of Garvey et al. However, in order to extract strangeness matrix elements of 
the nucleon from neutral current neutrino-nucleus scattering it is necessary to investigate 
the reliability of calculating the cross sections for this process by assuming finite strange 
quark form factors for the constituent nucleons. Previously, neutral and charged current 
neutrino-nucleus scattering have been investigated beyond the impulse approximation in 
but without strange quark form factors and MEG corrections. In addition, this and 
other proposals to extract strange quark form factors of the nucleon from neutrino- 

nucleus scattering involves kinematics ranging from low-energy inelastic scattering to the 
quasi-elastic region. Experience from electron scattering suggests that MEG corrections to 
neutrino-nucleus cross sections in this kinematic range might be important. 

The same LSND collaboration has recently announced their measurement of the cross 
section for the inclusive charged current reaction ^^C(z/^-, /i~)X near threshold |^ which is 
significantly lower than existing model predictions. This discrepancy is in some cases more 
than 100% ||3l| and it is conjectured that nuclear effects which are important in low-energy 
neutrino-nucleus interactions have been left out in these model calculations. One such 
effect is the MEG corrections to the impulse approximation and it is therefore interesting to 
explore whether two-body corrections to the one-body nucleon current will help to remedy 
this recently observed discrepancy. 

Previously, MEG corrections have been investigated only in neutrino-deuteron reac- 
tions using non-relativistic kinematics and without assuming any strange quark form factors 
32-p5|. In a recent letter [36| a method to construct relativistic MEG operators applicable 



to neutrino-nucleus scattering was presented taking into account the possible finite strange 
quark form factors of nucleons. This method treats all the SU{3) currents on the same 
footing and thus is able to estimate MEG corrections to electromagnetic as well as neutral 
and charged current processes simultaneously. It is also model independent in the sense that 
no underlying nucleon-nucleon interaction needs to be specified in order to construct the 
MEG operators. In this paper details of the method is presented together with examples 
of the use of the resulting MEG operators in electron and neutrino scatterings assuming 
nuclear matter and using the kinematics of the on-going LSND experiment. The focus of 
this paper is to investigate only the effects of exchange currents, and therefore a simple RFG 
model is used to model the target nucleus. In order to make a more realistic comparison 
with experiment additional nuclear effects must be incorporated in the description of the 
neutrino-nucleus scattering process. These effects will be considered in a forthcoming paper 

In the following section the problem is defined and empirical and theoretical motivations 
are presented for the use of the method originally developed by Ghemtob and Rho to con- 
struct MEG operators P3. Their formalism is then generalized to take into account finite 
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strange quark form factors and advantages of using the generalized method and assumptions 
made in constructing the MEC operators are discussed. Section III illustrates the usefulness 
of the generalized operators by evaluating two-body MEC corrections to quasi-elastic elec- 
tron scattering as well as neutral and charged current neutrino scattering reactions. Finally, 
the results are summarized in the concluding section accompanied by Appendicies which 
present all the necessary formalism needed to evaluate the nuclear response functions in 
neutrino-nucleus scattering together with some new technical details. 



II. EXCHANGE CURRENT OPERATORS FOR ALL CURRENTS 



A. Currents and Form Factors 



This section presents the construction of MEC operators applicable to both electron and 
neutrino-nucleus scattering by generalizing the method developed by Chemtob and Rho 
39| . As will be shown below, this generalized method is able to estimate exchange current 



corrections to any linear combination of SU (3) vector and axial currents given by 



a=0 



(2.1) 



where and A" are the SU{3) vector and axial vector currents of the nucleon, respectively, 
and k = is the four momentum of the incoming probe. The SU (3) singlet (a = 0) and 
octet (a = 1 — 8) currents are defined using the usual Pauli, (Q^), Dirac, F|'((5^), axial, 
G^(Q^), and induced axial H^^Q"^) on-s/ie// nucleon form factors where = —k"^ 



F;(fc)^. |[F°(g^)7, + F°(Q^)E, 



I 



A* 



I 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



Here i = 1 — > 8, A* are the usual Gell-Mann matrices normalized to Tr (A"A ) = 26ab, I is 
the identity matrix and the magnetic operator is defined as 



(2.6) 



with M being the free nucleon mass. Thus, the problem addressed in this paper is to 
estimate two-body MEC corrections to the general one-body SU (3) nucleon current shown 
in Eq. ( |2.1| ) when the nucleon is immersed in nuclear medium. Once this is accomplished, 
MEC operators corresponding to electromagnetic, weak axial, neutral and charged currents 
may be trivially constructed by taking the appropriate linear combinations. 

For example, the one-body nucleon electromagnetic current is recovered when = 1 and 
= 1/ \^ and by letting the remaining coefficients vanish, multiplied by the appropriate 
isoscalar and isovector electric charges, i.e. 
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(2.7) 



Similarly, one-body neutral, J^" and charged, J^*, 
linear combinations of vector and axial currents [|T9l 



currents are given by the following 



cos 6'c + 



4 ± 



sin 6*, 



(2.9) 



13 as well 



In these definitions of neutral and charged currents, 9w and 9c are the Weinberg and Cabbibo 
angles, respectively, and small QED, QCD and heavy quark corrections to 
as contributions from the charmed quarks to J"^^ 
Eq. 



■ ^ are ignored. Note that the third term in 
is proportional to the electromagnetic current, Eq. (^^71) , while the last two terms 



are referred to as the strange vector, V^, and axial, A* , currents of the nucleon. 



V V /3 
4^ = /lO 



(2.10) 
(2.11) 



For neutral current processes where only massless leptons are involved, the induced axial 
form factor does not contribute to the total cross section and therefore is not determined. 
At = 0, the strange quark Pauli form factor = — 2/\^Ff vanishes by definition 
and only F2 and are unknown among the form factors in Eqs. ( p.2|) -( p75|) . According 



to the Standard Model, these two form factors determine the strange quark magnetic, F2 = 
^2^(0), and axial, G% = G%{0), form factors of the nucleon. In this work all form factor 



parametrizations are taken from |^S[ which are the same ones used in Specifically, the 
dependence of -F|(Q^) and G^i^Q"^) are 



where 



(2.12) 
(2.13) 



F^m = 



F2°''(0) 



1M% 



)(1 



Ml 



(2.14) 



(2.15) 



In these definitions the vector and axial masses are set to My = 840 MeV and Ma = 1030 
MeV, respectively, and the octet form factors at = are -F|(0) = \/3/2{Kp + /t„) and 
G\{0) = VS/6{3F-D). 



5 



B. Soft Pion Exchange Dominance and the Chiral Filter Hypothesis 



Although empirical evidences abound suggesting that both one-body electromagnetic 
and weak axial currents are renormalized in nuclear medium by MEC, there still lacks a rig- 
orous theoretical understanding of the roles that different types of MEC might play in nuclei. 
For example, in electromagnetic processes any conserved transverse MEC consistent with a 
given N — N interaction is acceptable since it is not constrained by the Ward-Takahashi 
identity. A typical "brute force" approach of estimating MEC effects in electro-nuclear 
phenomena is to choose a model dependent N — N interaction in the one-boson exchange 
approximation and from this interaction construct the MEC operators with a longitudinal 
component satisfying the Ward-Takahashi identity and a corresponding conserved transverse 
component. In addition, in applications to many-body systems it is necessary to construct 
nuclear wave functions from the same N — N interaction including the Fock term to main- 
tain self-consistency. This latter requirement is often neglected without any justifications 
especially in relativistic calculations. The situation in the weak axial case is more ambiguous 
since there are no conservation laws to constrain the form of MEC operators except in the 
chiral limit. Thus, it is desirable to identify the important contributions from a multitude 
of two-body currents involving exchanges of different types of mesons by exploiting some 
underlying physical principles. Fortunately, there has been some promising experimental 
and theoretical progress over the past 20 years towards accomplishing this goal. 

In 1978, Kubodera, Delorme and Rho [0|, using the method developed by Chemtob and 



Rho to be described below, have argued that, in the absence of kinematical suppressions, 
MEC processes in low and intermediate energy nuclear phenomena are dominated by one- 
pion exchange whose production amplitude is evaluated in the soft-pion limit. In another 
words, they stressed that dominant contribution to the two-body correction arises from 
MEC which are consistent with nuclear force implied by chiral symmetry, other short-ranged 
meson exchanges being "filtered out" by the nuclear medium. Thus, the crucial element in 
their argument is the important role played by chiral symmetry in nuclei, a key symmetry 
manifest in QCD. Using this so-called chiral filter hypothesis, they arrived at the prediction 
of axial charge renormalization which subsequently was confirmed through the analysis of 
first forbidden /3-decay rates as mentioned in the introductory section. In general, Kubodera, 
Delorme and Rho found that MEC derived assuming soft-pion dominance approximation, 
hereafter referred to as soft-pion MEC, strongly renormalizes the time component of axial 
currents and the space component of electromagnetic currents, respectively. 

The idea of soft-pion dominance has been applied in the past to various low and in- 
termediate energy phenomena and proved to be a viable technique in estimating exchange 
current corrections. For example, in their analysis of first forbidden /5-decay rates Warbur- 
ton et al. compared the use of the soft-pion exchange dominance approximation to other 
MEC operators and found that both methods can reproduce the observed enhancement of 
the axial charge The successful application of soft-pion MEC operators here is not 

surprising since /3-decays typically involve small momentum transfers where the soft-pion 
dominance approximation is expected to apply. The real surprise came when the success 
of this approximation manifested itself in the electromagnetic sector. Initially, Riska and 
Brown used the soft-pion MEC operators as prescribed by Chemtob and Rho to explain 
the difference between the impulse approximation prediction and the measured threshold 
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radiative np capture rates. The same method was employed by Hockert et al. to calculate 
the cross section for the electrodisintegration of the deuteron which involved small energy 
{Enp ~ 3MeV) but large momentum transfers. They were able to reproduce the measured 
cross section up to momentum transfer of A;^ = 10 fm~^ using only the soft-pion MEC oper- 
ators. Addition of the A resonance contribution had little effect on their original correction 
to the impulse approximation. The same cross section was later measured at Saclay |^ 



extending the momentum transfer up to /c^ = 18 fm~ . The surprise came when the original 
prediction by Hockert et al. managed to reproduce the Saclay data up to fc^ = 15 fm~^ 
^ . In this case different corrections to the soft-pion MEC cancelled each other leaving the 
original two-body contribution to be the dominant correction to the impulse approximation. 
Thus, there are concrete empirical evidences to support the use of the soft-pion dominance 
approximation in low and intermediate energy electromagnetic and weak axial interactions 
in nuclei. 

Recently, Rho 02] has proposed an explanation for the success of the soft-pion MEC 



dominance based on Weinberg's derivation of nuclear forces from chiral Lagrangians |43,44 



Using chiral power counting arguments Rho has shown that to the leading order, i.e. at 
the tree level, the short range part of two-body MEC corresponding to a nuclear force 
predicted by a given chiral Lagrangian is considerably suppressed. Thus, the dominant 
contribution to the two-body correction to the impulse approximation comes from the long 
ranged part of the exchange current represented by the soft-pion exchange. Subsequently, 
Park, Towner and Kubodera have calculated next-to-leading order corrections to the 
axial charge MEC operators beyond the soft-pion dominance approximation using heavy 
fermion chiral perturbation theory. They found that loop corrections to the soft-pion MEC 
operators are of the order of 10%, and argued that their results are consistent with the 
claims of Warburton et al. and support the chiral filtering conjecture. Thus, not only are 
there empirical evidences suggesting the renormalization of electromagnetic and weak axial 
currents by MEC in nuclei, but there also exists theoretical motivations to believe that this 
renormalization is dominated by soft-pion exchange interaction, at least up to momentum 
transfers of about one GeV El|,E^. 



C. Soft Pion MEC Operators 

The method of Chemtob and Rho to construct soft-pion MEC operators is based 
on soft-pion theorems and current algebra techniques pioneered by Adler Here, this 
method is generalized to accommodate all the SU{3) currents appearing in Eq. ( |2.1| ) and the 
main advantage of using this technique to estimate MEC corrections to neutrino-nucleus 
scattering is pointed out. Since the derivation of soft-pion MEC operators may be found in 



the original works of Adler [46| and Chemtob and Rho |39], and the generalization to SU{3) 
being straightforward, most of the technical details are relegated to Appendix B. 

In the soft-pion dominance approximation, the operator representing an exchange of a 
pion between two nucleons is written as products of the pion production amplitude by an 
external current off the first nucleon, the pion propagator and the matrix element for pion 
absorption by the second nucleon 

•^u(^; PiX^ PiXi Pf,i] Pf,2)ex = 



7 



(2vr: 



-6\Pj,, + P,,2 + k - Pf,i - PF,2){N{PF^)7r\q)\j;{k)\N{Pj^,)) 



q2 _ ^2 



{N{Pp,2)\Jt{q)\N{Pi,2)) + il 



(2.16) 



Here the matrix element (A^(PF^i)7r''(g)| J^(/c)|A^(P/ i)) is the amphtude for pion production 
off a nucleon by an external vector or axial SU{3) current J^{k), J^{k) + iV(Pf^i) — »• 7r^(g) + 
N{Pf^i), where q = is the four momenta of the produced pion. a and b are SU{3) indices 
with a = — ^ 8 to accommodate all the SU{3) currents and 6 = 1, 2 or 3 for pion production. 
Similarly, {N{PF^2)\J^{q)\N{Pi^2)) is the matrix element for pion absorption by a nucleon 
using the pseudoscalar ttA^A^ coupling 



{NiPF^2)\Jtiq)\NiPl,2)) = g.NN (iV(P^,2)|75AliV(P,,2)) 



(2.17) 



Note the absence ofirNN form factor in Eg. /\2.17{ ). It is quite remarkable that the present 
method can describe the electrodisintegration data involving momentum transfers of up to 
about k"^ = 15 fm without the use of any tcNN form factors in the tcNN absorption vertex. 

To construct soft-pion MEC operators it is necessary to know the pion production am- 
plitude in Eq. ( ^.16 ) in the soft-pion limit of g 0. In this limiting procedure it is necessary 
to first take the spatial part of the four vector to zero in order to select the long range, i e. 
the S-wave, part of the N — N interaction and then take the chiral limit of go ~^ 0. The 



resulting amplitude, originally derived by Adler using the PCAC hypothesis 



= m^P^vr" 



(2.18) 



and used by Chemtob and Rho in has the following form when generalized to the SU (3) 
formalism 



lim{N{PF)n\q)\j;{k)\N{Pj)) = —J d^x lim(-0(iV(P^)|T (^(a;) J,"(0)) |iV(P,)) (2.19) 



(iV(P^)l Ql{x),rJO) _\N{Pj)) 



xo=0 



Here Q'^{x) = J (I^xAqIx) is the axial charge. As shown in the only contributions to 
the first term in the soft-pion limit come from pole diagrams where the axial current is 
inserted in the external lines in the amplitude (A^(Pp)| J°(/c)|A^(P/)) and thus behaving as 
1/g^. The second term may be simplified by using the SU{3) ® SU{3) current algebra 



g«(x),4(o)]^^^ = z/,,eV;^(o) 



(2.20) 
(2.21) 



and has no contributions from singlet currents unlike in the first term where both SU (3) 
singlet and octet can contribute. Therefore, the amplitude for soft-pion emission in the re- 
action N{Pj) — > N{Pp), in the presence of perturbation J^{k), may be expressed in terms of 

two matrix elements {N{PF)\J?,{k)\N{Pi)) and {N{Pf)\ \qI{x) , J",{0)\ \N{Pj)). Since 

J" may be any one of vector or axial SU{3) currents, Eq. ( p.l6|) may simultaneously be 
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applied to all the components of the general one-body SU{3) current in Eq. ( |2.1|) , and 
specifically to the electromagnetic current of Eq. ( p.7|) as well as to neutral and charged 
currents of Eqs. ( |2.8| ) and ( |2.9| ), respectively. Thus, it is the use of SU{3) ® SU{3) cur- 
rent algebra in Eq. (|2.2CI| ), which rotates around the vector and axial currents, that makes 
the generalized method of Chemtob and Rho particularly suitable to approximate MEC 
corrections in neutrino-nucleus scattering at low and intermediate energies. In addition, 
because the soft-pion limit is taken no N — N interaction needs to be specified and the 
present method of constructing MEC operators only requires the currents in the impulse 
approximation as inputs. In this sense the present approach to constructing MEC operators 
is model independent. Note also that the present method is valid to all orders in Qt^nn 
in one of the irNN verticies since the pion production amplitude Eq. ( p.20|) is evaluated 
non-perturbatively [^j .f] 

It is useful to discuss the approximations made in calculating the soft-pion production 
amplitudes. Since the original application of soft-pion theorems has been on pion photo- 
production processes, these amplitudes have been evaluated by assuming that the initial 
and final nucleons are on their mass shell. However, in order to construct MEC operators 
in nuclei consistently it is necessary to take both the initial and the final nucleons off their 
mass shell. This involves density dependent off-shell parametrizations of nucleon currents 
which are not known. In fact, a fully consistent many-body description of in-medium 
nucleon electromagnetic or weak form factors has never been presented. Considering this 
lack of understanding of off-shell modification of nucleon form factors and currents in nuclear 
medium, the most reasonable approximation to make is the use of on-shell kinematics and 
parametrizations of nucleon form factors as has been done in all previous works involving 
MEC in nuclei. This implies that the nucleons are assumed to obey the free Dirac equation, 
and therefore in the derivation of the pion production amplitude "^u^p) has been replaced 
with Mu{p), where u{p) is the nucleon spinor with three momentum p and mass M. This 
on-shell approximation is used from the very beginning, even before considering the soft- 
pion limit, and is consistent with using the free RFC model of the nucleus to calculate the 
cross sections where the constituent nucleons are assumed to be on-shell. 

Related to the derivation of soft-pion production amplitudes, Eq. (|2.20| ), is the use of 
PCAC and pion off-shell extrapolation. When using the PCAC relation, it is assumed that 
the matrix element of the divergence of the axial current varies smoothly and slowly with 
the pion mass q^. In addition it is also assumed that this matrix element is approximately 
proportional to the pion field and that any dependence on higher order non-linear pion field 
terms are negligible. This assumption allows one to make a connection between the results 
obtained with g ^ with a more realistic value of q^. Since a zero four momentum pion 
is not a physical object, this slow variation assumption is always needed to compare the 
soft-pion predictions with experiment. Finally, the soft-pion production amplitudes are 
valid only to zeroth order in q. This means that if there are any processes contributing to 



^The expressions for the soft-pion production amplitudes presented in Appendix B are propor- 
tional to QttNN- However, this does not mean that the amplitudes are evaluated to first order in 
9nNN- The ttNN coupling constant was introduced through the use of the Goldberger-Trieman 
relation |46|. 
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the pion production which, for kinematical reasons, are of first order in g, the apphcations of 
current algebra and soft-pion techniques used here give no information about them. In order 
to describe these processes, it is necessary to make a model dependent analysis of specific 
reactions which requires the introduction of a — interaction. 

Relativistic SU (3) soft-pion MEC operators are constructed in analogy with the method 
outlined in and are listed in Appendix B. The conservation of the vector current has 



been checked analytically using the prescription discussed in the original paper by Adler 
and verified numerically . To illustrate the usefulness of the present technique, soft-pion 
MEC corrections are applied to quasi-elastic electromagnetic, neutral and charged current 
interactions simultaneously in the following section. In this paper a simple RFG model 
formalism ||48|j49(] is used for all calculations without binding energy corrections {B = 0). 
Finite nucleus effects, final state interactions and other density dependent nuclear medium 
effects are not considered on purpose in order to clearly isolate the effects of soft-pion MEC 
in many body systems. 



III. RESULTS 

A. Quasi Elastic Electron Scattering 

As mentioned in Section IIA, the third term in the expression for the one-body neutral 
current of the nucleon, Eq. (|2.8|), is proportional to the electromagnetic current, (xV^ + 
-^V^^. Thus, estimates of MEC effects in electron scattering can automatically be extracted 
when calculating two-body corrections to neutral current processes using the formalism 
under consideration. The successful application of the soft-pion MEC approximation to 
the electrodisintegration of the deuteron has been described in the previous section. Here, 
the same method is applied to quasi-elastic electron scattering off heavier nuclei where 
many-body effects, not present in the reaction involving the deuteron, are expected to play 
important roles. 

In Figures la and lb, separated longitudinal, Rl{^^, \k |), and transverse, Rt{uj, \k |), 
response functions for a typical quasi-elastic inclusive ^^C(e,e') reaction are shown as a 
function of energy transfer u both in the impulse approximation and with soft-MEC cor- 

— * 

rections. A Fermi momentum of kp = 225 MeV and a fixed three momentum transfer of |fc | 
= 400 MeV are used in the RFG model of the nucleus without any binding energy correc- 
tions {B = 0) so that the nucleons in the target nucleus are on-shell. Since the momentum 
transfer is almost twice the Fermi momentum, the influence of Pauli blocking should be 
small in this kinematical range and indeed, as shown in the figures, the small effect of Pauli 



blocking is manifested in the linear dependence of the response functions for small uj |^ 
Thus, the RFG model of the nucleus is adequate as a first approximation to test the soft- 
pion dominance approximation. Only lp\h final states are considered when evaluating the 
MEC matrix elements since they are the dominant two-body contributions to the response 
functions in the quasi-elastic region. A more quantitative two-body MEC corrections to 
electron scattering at these energies, and especially in the dip region, would require model 
dependent A propagation and pion production processes as pointed out in ||9|JT0|1. Also 



shown in the figures for reference are the experimental data for the same inclusive electron 
scattering reaction measured at Saclay However, no attempt has been made to fit the 
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data by varying kp and B, since the vector currents will no longer be conserved with a finite 
binding energy correction and the target nucleons will be off their mass shell. 

As shown in Figure la, the MEC correction to the impulse approximation in the longi- 
tudinal response results in a small reduction of the quasi-elastic peak with no appreciable 
change in the peak position. However, this reduction of the impulse approximation result is 
too small to describe the magnitude of the quasi-elastic peak observed in the Saclay data. 
This feature persists even if kp and B are varied in an attempt for a better fit to the data. 
The failure of the free Fermi Gas model to describe the observed longitudinal response in 
quasi-elastic inclusive electron scattering for a wide range of nuclei is well-known |]^, and 
additional nuclear effects must be incorporated in order to improve the present model cal- 
culation. Also shown in the figure is the individual soft-pion MEC contribution which is 
small and negative resulting in the small reduction of the quasi-elastic peak. 

For the transverse response shown in Figure lb, the soft-pion MEC correction can in- 
crease the magnitude of the quasi-elastic peak by about 20% relative to the impulse approx- 
imation results and shift the peak position to a lower value of energy transfer u by about 
20 MeV. These effects of the two-body current correction may be understood by examining 
the individual soft-pion MEC contribution which causes both constructive and destructive 
interferences depending on the value of the energy transfer. This contribution is positive 
between a; = and 125 MeV and negative thereafter and results in an increase of the mag- 
nitude of the quasi-elastic peak and a shift of the peak position to a lower value of uj. The 
qualitative agreement with data is better than in the longitudinal response although the 
observed peak position can not be reproduced with zero binding energy correction {B = 0). 
A finite value of the binding energy will help to shift the quasi-elastic peak towards the 
observed position but, as mentioned above, the total electromagnetic current will no longer 
be conserved in this case. Thus, two-body soft-pion MEC corrections to inclusive electron 
scattering response functions tend to slightly soften the longitudinal response and leads to 
about 20% increase in the transverse response in the free RFC model. 

It is interesting to compare the present results with numerous model calculations of Rl 
and Rt in the literature. However, to make any quantitative comparisons with other works 
are very difficult since many model dependent assumptions have been made for each of 
the calculations which are hard to disentangle. For example, Kohno and Ohtsuka also 
investigate the MEC corrections to the inclusive electron scattering response functions in 
the quasi-elastic region considering only the Iplh final states. In contrast to the present 
work, they work in the non-relativistic approximation and ignore contributions of order 
greater than glj^j^. The method adopted in this paper is fully relativistic and one of the 
ttNN verticies is valid to all orders in Qj^nn- Also, they include the A excitation current 
which has a strong model dependence as well as the "pion-in-fiight" MEC, referred to as 
"pionic" MEC in their paper. As discussed in Appendix B, the pion-in-fiight diagram does 
not contribute to the (e, e') response functions since it is proportional to the four momentum 
transfer of the leptonic probe, fc^, in the soft-pion limit. Nevertheless, there is qualitative 
agreement between the soft-pion MEC contribution shown in Figures la and lb and their 



corresponding "pion pair" MEC results presented in Figure 2b of [^. It is interesting to 



note that in their calculation for the longitudinal response, there is a cancellation between 
the "pionic" and A excitation MEC corrections leaving the "pair" current the dominant 
correction to the impulse approximation as in Figure la of the present paper. This is not 
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true for their transverse response where Kohno and Ohtsuka find that the total effect of the 
two-body corrections is to decrease the impulse approximation result which is the opposite 
effect shown in Figure lb. The reason for this discrepancy is due to the large negative 
contribution from their A excitation current which is not included in the present work. It 
should be noted that another RFG model calculation of Rt that includes the same types 
of MEC corrections as in Kohno and Ohtsuka [0 , but considers only the 2p2h final states, 
obtains an increase in the transverse response since the 2p2h matrix elements are added 
incoherently to the impulse approximation result. Thus, unlike the case of the deuteron, 
it is difficult to make a definitive statement on the success of the application of soft-pion 
MEC dominance approximation in electron scattering off nuclei by comparing with other 
works, except to say that the results presented above are not inconsistent with other model 
calculations. 



B. Quasi— Elastic Charged Current Neutrino Scattering 



The one-body charged current, Eq. ( p^.9D , is the next simplest linear combination of 
SU{3) vector and axial currents of physical interest after the electromagnetic current dis- 
cussed above. Since the nucleon has no net strangeness {Ff{Q^ = 0) =0), only the term pro- 
portional to cos6'c in Eq. ( pl9| ) contributes to the charged current neutrino-nucleus scattering 
cross sections. Although no direct information regarding the strange quark form factors of 
the nucleon may be obtained in charged current processes, these reactions are nevertheless 
important in extracting the value of the axial mass parameter Ma which appears in the 
parametrization of axial form factors. Recently, the LSND collaboration has reported on 
their measurement of neutrino flux-averaged inclusive ^"^Clu^- , fi~)X cross section of (8.3 
± 0.7 stat. ± 1.6 syst.)xlO~^° cm^ in the neutrino energy region of 123.7 < < 280 MeV 



with a flux-weighted average oi < Ey > = 180 MeV |3^. This value is substantially smaller 
than an earlier measurement, which used a different neutrino energy distribution, of (15.9 
± 2.6 stat. ± 3.7 syst.)xlO~^^ cm^ [^| as well as predictions from model calculations ||31|| . 



Because these cross sections are measured just above the muon threshold energy, it is ex- 
pected that nuclear effects will substantially modify the impulse approximation predictions. 
MEC corrections is one such effect and its impact on the impulse approximation calculation 
is examined in this subsection using the soft-pion dominance formalism. 

In contrast to the charged current reaction involving the muon-neutrino, it is interesting 
to note that there is good agreement between model predictions and the measured cross sec- 
tion for the exclusive charged current reaction ^'^C{i'e,e~)Ng s. involving e/ectron-neutrinos 
53. From the computational point of view, the main difference between electron- and 



muon-neutrino induced charged current reactions is that in the former case, the mass of 
the outgoing electron is negligible while the muon mass has to be explicitly included for the 
latter case. This means that in the charged current reaction ^^C(z/^- , yU~)X, the term pro- 
portional to the induced axial form factor, if^(/c^), will contribute to the total cross section. 
In the present formalism, the induced axial form factor is given by PCAC in the pion pole 
dominance 

HAiP) = Ga{Q)t^^ (3.1) 
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where Ga(0) ^ Qa/'^' = + D) is the axial couphng constant. However, because of the 

presence of the HAik"^) term there is an ambiguity when applying current conservation to 
the vector currents. The prescription given by Adler in has been used here to insure 
vector current conservation and details are given in Appendix B. 

Figure 2 shows the total cross section for the inclusive charged current process 
^^C(z/^- , for several nuclear densities obtained by folding the LSND neutrino energy 
distribution [^U|,^. As in the electron scattering case, the RFG model without binding en- 
ergy corrections is used to model the nucleus and only Iplh final states are considered when 
taking matrix elements of two-body operators since the phase space for 2p2h final states 
should be suppressed due to the rather low energy neutrino beam [Q. However, because the 
flux-weighted average of the neutrino beam energy is < -E^, > = 180 MeV, the three momen- 
tum transfers are twice smaller than in the electron scattering reaction discussed above and 
consequently, the effect of Pauli blocking becomes more important in this charged current 
reaction. In the RFG model of the nucleus an increase in nuclear density is equivalent to 
an increase in the Fermi momentum. Thus, with a constant neutrino beam energy the net 
effect of Pauli blocking is the decrease of the total cross-section per nucleon with increasing 
nuclear density as shown in Figure 2. In the impulse approximation, the total cross section 
decreases from 24 to 13 (xlO~^°cm^) as the Fermi momentum is varied from kp = 220 to 
300 MeV. As shown in the figure, the inclusion of soft-pion MEG corrections reduces the 
impulse approximation results by 5 to 10% as the Fermi momentum is increased from 220 to 
300 MeV. For kp = 225 MeV, which is the usual value used for ^^C, the total cross section is 
reduced from 24.1 to 22.7 (xlO~^''cm^). This reduction is not enough to explain the recently 
measured value reported by the LSND collaboration. 

In Figure 3a, the ^^C(z/^-, cross section is shown as a function of neutrino energy 
for kp = 225 MeV including Goulomb corrections for the outgoing muon. Gorrections from 
two-body currents to the impulse approximation are very small and difficult to observe. 
Note that the muon production threshold is found to be around 107 MeV compared to 
the experimental threshold of 124 MeV. This is because the present RFG calculation does 
not include any binding energy corrections which will ruin the vector current conservation. 
Gorresponding calculation using non-relativistic continuum RPA formalism obtained a 
smaller cross section which increases more slowly with the neutrino energy. The effect of 
soft-pion MEG is more evident when the differential ^^C(z/^-, /i~)X cross section folded with 
the experimental neutrino energy distribution is plotted against the muon kinetic energy, i?^, 
as shown in Figure 3b. Gorrections from two-body currents are largest near the cross section 
peak around = 23 MeV where the impulse approximation result of 43.7 x 10~^^ cm^/MeV 
is reduced to 40.1 x 10~^^ cm^/MeV. However, these corrections become less important with 
increasing muon kinetic energy and for E^^ > 60 MeV there is hardly any change from the 
impulse approximation result. The results shown in Figure 3b are qualitatively consistent 
with non-relativistic Fermi Gas and continuum RPA results shown in |^ although various 
model dependent assumptions, such as a different neutrino energy distribution, prevents from 
a direct quantitative comparison between these results. Thus, although the two-body soft- 
pion MEG corrections help to reduce the impulse approximation prediction of the exclusive 
charged current reaction ^^C(t'^-, towards the observed value, this reduction is not 
enough and additional nuclear structure effects are required to further lower the total cross 
section for this exclusive process ||5^J57[] . These additional effects will be incorporated in a 
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future work 



C. Quasi Elastic Neutral Current Neutrino Scattering 

The final application of the soft-pion MEC correction is in the quasi-elastic neutral 
current neutrino-nucleus scattering where the one-body neutral current of the constituent 
nucleon is given by Eq. ( |2.8| ). This particular linear combination of the general SU (3) current 
of Eq. is interesting since, unlike the electromagnetic and charged currents, it includes 
the SU{3) vector and axial singlet currents leading to a possible strange quark contribution to 
the neutrino-nucleus scattering reaction. Also, in contrast to the charged current reactions, 
the term proportional to the induced axial form factor Ha does not contribute to the neutral 
current neutrino-nucleus cross section since the leptons involved in the scattering process 
are massless. In Eq. (|2.8| ), the strange quark electric form factor is assumed to be 
vanishingly small for small values of virtuality that are of interest here, and therefore only 
magnetic, F2 and axial G\ strange quark form factors of the nucleon are used to describe the 
contributions from the strange sea quarks in neutral current neutrino-nucleus scattering. In 
the following discussion these two strange quark form factors are treated as input parameters 
and no attempt has been made to determine their values from model calculations. As in the 
preceeding discussions on electromagnetic and charged current reactions, all the calculations 
are performed using the RFG model without binding energy corrections. 

In Figure 4, differential cross sections for the ^^C(z/, i/'p) reaction, da/dEp, where Ep 
is the total energy of the ejectile, are shown in the quasi-elastic region as functions of the 
ejected proton's kinetic energy, Tp, for several values of Fermi momenta. The incident 
neutrino energy is assumed to be 200 MeV and the values used for the strange magnetic and 
axial quark form factors are F2 = —0.21 and G\ = —0.19, respectively. The long dashed 
curves are the impulse approximation results while the solid curves have been obtained with 
the soft-pion MEC corrections, and as in the charged current case, only Iplh final states 
have been taken into account. In the impulse approximation, the magnitude of the peak 
decreases from about 70 to 53 (xlO~^^ cm^/MeV) as kp is varied from 200 to 350 MeV. 
However, this decrease is accompanied by a redistribution of the strength to higher values of 
Tp in such a way that the area under the differential cross section remains approximately a 
constant since it is roughly proportional to the number of nucleons which is fixed. Note that 
because the incident neutrino energy is constant and small, the position of the maximum of 
the differential cross section is approximately the same at around Tp=30 MeV despite the 
decrease in the peak magnitude. This means that the effect of Pauli blocking becomes very 
important when computing the total cross section obtained by integrating the differential 
cross section da/dEp for Ep > Epermi, thereby cutting off the contributions around the 
peak magnitude. This total cross section decreases noticeably with increasing kp, as it was 
the case for the total charged current cross section shown in Figure 2. 

It is evident from the figures that the soft-pion MEC effects are sensitive to the nuclear 
density in this model calculation. For kp around 200 MeV, there is very little difference 
between the impulse approximation and the MEC corrected results due to an important 
cancellation to be discussed below. This cancellation becomes less complete with increasing 
kp and the two-body current effects become stronger and decreases the impulse approxi- 
mation results with increasing nuclear density as shown in the figures. This effect of the 
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soft-pion MEC on the impulse approximation persists when the strange quark form factors 
are assumed to vanish {F2 = G% = 0) in the model calculation. However, without any 
strange quark form factors there is a large reduction in the differential cross section for pro- 
tons as shown in Figure 5. In this figure kp is set to 225 MeV and there is little correction 
coming from the exchange currents in the soft-pion approximation for this value of Fermi 
momentum. Note that the magnitude of the quasi-elastic peak is reduced from about 70 
to 50 xlO~^^ cm^/MeV when the strange quark magnetic and axial form factors are varied 
from F2 = —0.21 to F| = and G\ = —0.19 to G\ = 0, respectively. For the neutral 
current induced neutron knockout reaction, ^^C(z/, i/'ra), a finite strange quark form factor 
decreases the differential cross section relative to that obtained without any strange quark 



distribution for nucleons 22,25 



It is instructive to understand why the soft-pion exchange current corrections are small 
for neutral current scattering for nuclear densities corresponding to kp ~ 200 MeV. As 
shown in Appendix A Eq. ( |A38| ), the differential cross section per nucleon for neutral current 
scattering in the RFG model is given by a linear combination of three structure functions 
Wl, Wt and Wtt as 

'd(7^^^\ ( Svr \ (7^^ r 1 

defd{cos9) ojlWl + ojtWt + ojtt'Wtt' (3.2) 



^^^/rfgm \k 



Here a;^, ujt and ujtt' are the kinematical coefficients corresponding to the longitudinal Wl, 
transverse Wt and transverse-transverse Wtt' structure functions, respectively. In Figure 6 
the individual contributions from these three structure functions to the differential cross 
section are shown using the same parameters as in Figure 5. Note that all three contri- 
butions to da/dEp are positive, the dominant one coming from the transverse term ujtWt 
followed by the transverse-transverse utt'Wtt' and the longitudinal ujlWl terms, the latter 
contributing negligiblely to the differential cross section. Note that the soft-pion MEC cor- 
rection increases the ujtt'Wtt' term and decreases the uopWp term relative to their impulse 
approximation results leading to a cancellation of MEC effects for this particular value of 
Fermi momentum. This cancellation between the transverse and transverse-transverse con- 
tributions becomes less and less complete as the nuclear density is increased and the effect of 
two-body correction becomes bigger for higher densities as shown in Figure 4. Eq. ( |3.2|) can 
also be used to understand the behaviour of the differential cross section for anti-neutrino 
scattering relative to the neutrino induced nucleon knockout rates as follows. 

Figure 7 shows the da/dEp for the neutral current anti-neutrino scattering reaction 
^^C(z/,z/'p) and the corresponding neutrino reaction both of which are obtained again by 
using the same input parameters as in Figure 5. In the impulse approximation, the mag- 
nitude of the quasi-elastic peak for anti-neutrino scattering is approximately half that 
of the neutrino scattering, the former value being 38 xlO~^^ cm^/MeV compared to 70 
xlO^^^ cm^/MeV for the latter. Also, the position of the peak is shifted towards a lower 
value of Tp by about 10 MeV relative to the neutrino scattering case. The effect of soft-pion 
MEC corrections to the impulse approximation result for anti-neutrino scattering is notice- 
ably larger than the corresponding correction for neutrino scattering for nuclear densities 
corresponding to kp = 225 MeV as shown in the figure. For the ^^C(i^, z/'p) reaction, the im- 
pulse approximation result for the magnitude of the quasi-elastic peak is reduced from 38 to 
30 xlO~^^ cm^/MeV, while there is no such noticeable reduction for the neutrino scattering. 
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Thus, not only the anti-neutrino scattering differential cross section greatly reduced relative 
to neutrino scattering, but the effect of soft-pion MEC on the impulse approximation result 
is much bigger reducing further the da/dEp by about 15%. 

This difference in the neutrino and anti-neutrino differential cross sections may be under- 
stood by once again examining Figure 6 and noting the definitions of the structure functions 
Wi and their corresponding kinematical coefficients Ui given in Eqs. ( |A15|) to ( |A29|) of Ap- 
pendix A. The only difference between the longitudinal, transverse and transverse-transverse 
structure functions and kinematical coefficients entering in Eq. (|3.2|) for the neutrino and 
anti-neutrino scattering is in the sign of ujtt' given by Eq. ( |A28 ). Therefore, the graph corre- 
sponding to Figure 6 for anti-neutrino scattering will have a negative transverse-transverse 
contribution which partially cancels the dominant transverse contribution yielding the small 
differential cross section as shown in Figure 7. In addition, because the ujtt'Wtt' term is 
now negative, the MEC correction to this term will further decrease the correction to the 
transverse term resulting in a non-negligible two-body correction to the differential cross 
section not observed in the neutrino scattering case. Hence although the differential cross 
section for the anti-neutrino reaction ^^C(z/, i^'p) is about half that of the corresponding 
neutrino scattering reaction, its impulse approximation result is subject to non-negligible 
soft-pion MEC corrections. 

In the Introduction it was mentioned that the ratio of proton-to-neutron yields in quasi- 
elastic neutral current scattering, R{p/n), is currently being measured in the LSND exper- 
iment to extract the value of the strange quark axial form factor, G%. In this experiment 
the probing neutrinos and anti-neutrinos are incident on a tank of mineral water composed 



of hydrogen and carbon molecules |]22|. The amount of kinetic energy carried by the struck 
nucleon is measured by the detectors surrounding the tank of mineral water and R{p/n) 
is determined as a function of the detected ejectile's kinetic energy, Tp. The proton-to- 
neutron ratio is then integrated over Tp and compared with model predictions which gives 
the integrated R{p/n) as a function of G% for several values of the strange quark magnetic 
form factor as shown in Figure 8. However, in practice there is an experimental cutoff in the 
range of Tp, given by 50 MeV < Tp < 120 MeV, in order to make sure that the detected 
protons are knocked out from the carbon nucleus. The RFC model calculation of the ratio 
of integrated proton-to-neutron yields of Figure 8 is qualitatively similar to the ratios pre- 
dicted by a non-relativistic RPA calculation |^ and shows only about a few percent change 
from the impulse approximation results in the ratio for the neutrino scattering while this 
change is about 10% for anti-neutrinos. The latter change originates in the non-negligible 
two-body correction to the impulse approximation in the anti-neutrino scattering reaction 
explained in the preceeding paragraph. In both cases, the two-body corrections tend to 
decrease the integrated ratio and this effect increases as G\ becomes less and less negative 
for a given value of For example, if G% ~ —0.1, the soft-pion MEC corrections are neg- 
ligible for the neutrino scattering R{p/n) for all physically relevant values of F2 but must 
be taken into account in analyzing the anti-neutrino scattering ratio. Finally, final state 
interaction effects should be examined within a relativistic framework, although it is known 
from a non-relativistic calculation that the effects on the integrated ratios are almost 
negligible. 
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IV. DISCUSSION 



In this work relativistic two-body MEC corrections have been apphed simultaneously 
to quasi-elastic electromagnetic and neutrino scattering reactions using a single unified 
formalism obtained by generalizing the method developed by Chemtob and Rho [^. The 
basic strategy employed here is to use the chiral filtering hypothesis which assumes 



that the two-body MEC correction in nuclei is dominated by the exchange of a single pion 
whose production amplitude is evaluated in the soft-pion limit. By assuming PCAC and 
using the SU{3) ® SU{3) current algebra this soft-pion dominance approximation allows 
one to express the pion production amplitudes in terms of SU (3) vector and axial currents 
interacting with the nucleon. Thus the resulting MEC operators are expressed in terms 
of currents used in the impulse approximation without any references to model dependent 
N — N interaction. These soft-pion MEC operators arise from the longest range {i.e. the 
S-wave) part of the N — N interaction and PCAC is used to relate the soft-pion predictions 
to experiment. This is the reason why it is necessary to take the spatial part of the pion 
four momentum to zero first when taking the limit. The generalized version of 

the method of Chemtob and Rho used in this work is so far the most economical way 
to estimate exchange current corrections to low and intermediate energy neutrino-nucleus 
scattering since it treats all the SU{3) vector and axial currents in Eq. ( |2.1| ) on the same 
footing, and is a natural way to accommodate the strange quark vector and axial currents 
of the nucleon allowed by the Standard Model. 

Previously, the non-relativistic and SU{2) version of the soft-pion MEC operators have 
been successfully applied to first forbidden /3-decays |ll3| , [l5H r7| and to electromagnetic and 
neutrino scattering off few nucleon systems [^,0,^,^ including processes involving large 
momentum transfers 0. This last phenomenological success came as a real surprise since 
the soft-pion MEC operators were expected to be applicable only to those processes in- 
volving small momentum transfers such as in /3-decays. It led to speculations that the 
physics in nuclei is dominated by processes dictated by chiral symmetry and a theoretical 
justification was proposed to explain the phenomenological success by extending the 
work of Weinberg |Q on nuclear forces derived from chiral Lagrangians [^]. Motivated by 
these developments, fully relativistic SU{3) soft-pion MEC operators have been used to es- 
timate two-body corrections to the impulse approximation in quasi-elastic electromagnetic 
and neutrino scattering reactions assuming a finite magnetic and axial strange quark form 
factors of the nucleon. In order to clearly isolated the effect of two-body corrections in many 
body systems, a simple free RFC model of the target nucleus was used without any other 
medium dependent corrections. 

The inclusive (e, e') longitudinal and transverse response functions in the quasi-elastic 
region obtained using MEC corrections show a small reduction of the impulse approximation 
result for the longitudinal channel while a substantial increase was observed for the trans- 
verse response. The former result is in qualitative agreement with a previous work using 
non-relativistic kinematics [E^i, while the relativistic MEC corrections tend to increase the 



response function in the transverse channel in contradiction to the results of but in 
agreement with another relativistic calculation It should be stressed that in almost 

all calculations of MEC corrections in electron scattering reactions, such as in |I0|J5^ , the 
one-pion exchange operators are of order g^j^^^ whereas the MEC operators used in this 
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work is valid to all orders in Qt^nn at one of the irNN verticies since the soft-pion produc- 



tion amplitude is evaluated non-perturbatively [^. Thus it is difficult, if not impossible 



to make any quantitative comparisons with other works dealing with MEC corrections in 
electron scattering especially if additional medium corrections are present in a given calcula- 
tion. However, in order to make a more quantitative prediction it is necessary to introduce 
these nuclear medium effects, such as density dependent nucleon mass and form factors, 
as well as MEC involving the excitation of the A resonance. For these reasons, the most 
successful application of soft-pion MEC correction in electromagnetic interactions is still in 
the electrodisintegration of the deuteron investigated over 20 years ago. 

The charged current neutrino-nucleus reaction investigated in the present work is the 
inclusive ^^C(z/^- , reaction where the total cross section has recently been measured 



by the LSND experiment I^Oj. In the RFC model calculation, the two-body soft-pion 
exchange correction reduces the impulse approximation result for the cross section by 5 
to 10 % as the Fermi momentum is increased from 200 to 300 MeV. For kp = 225 MeV, 
a typical Fermi momentum used for ^^C, the total cross section is reduced from 24.1 to 
22.7 (xlO~^'^cm^) which is not enough to explain the recently measured value reported by 



the LSND collaboration of (8.3 ± 0.7 stat. ± 1.6 syst.) x lO'^^cm^ g. Recall that the 



previously reported measurement is (15.9 ± 2.6 stat. ± 3.7 syst.) x 10~^^cm^ [B3 which was 
measured using a different neutrino energy distribution than in the LSND experiment. If the 
LSND measurement is correct, then some important density dependent physics is responsible 
for the large reduction of the impulse approximation result, a situation reminiscent of the 
"missing strength" problem in the longitudinal (e, e') response functions of the last decade 

The quasi-elastic neutral current induced nucleon knockout reaction ^^C(z/, u'N) is cur- 
rently being investigated by the LSND collaboration in an effort to extract by measuring 
the ratio of proton-to-neutron yields. For the proton knockout reaction, the impulse ap- 
proximation result for the differential cross section is found to be sensitive to the values of 
strange quark form factors of the nucleon, while the two-body soft-pion correction becomes 
more important as the nuclear density is increased. However, for densities appropriate for 
the MEC corrections to the structure functions cancel each other resulting in small 
changes from the impulse approximation results. Because of this cancellation the ratio of 
proton-to-neutron yield for neutrino scattering does not suffer from two-body corrections 
and, at the present level of sophistication, the impulse approximation result is sufficient to 
describe the ratio R{p/n). The corresponding knockout reaction involving anti-neutrinos 
exhibits a much smaller differential cross section than in the neutrino scattering case and a 
non-negligible MEC correction to the impulse approximation result. These differences may 
be understood by once again examining the contributions from different structure functions 
to the differential cross section. In this case, one of the kinematical coefficients correspond- 
ing to a structure function changes sign relative to the neutrino scattering case. As a result, 
there is a cancellation in the impulse approximation results leading to a reduced differential 
cross section but an enhancement in the magnitude of the MEC correction which amounts to 
about 10-15% of the impulse approximation result depending on the value of G% However, 
in order to make a concrete statement about the extraction of the strange quark axial form 
factor from the measured R{p/n), additional nuclear effects, such as final state interactions, 
need to be considered in addition to MEC corrections. 
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In this paper, motivated by the LSND neutrino scattering experiment, the quasi-elastic 
region in electron and neutrino scattering was chosen to apply the soft-pion MEC correc- 
tions. This kinematic regime probably represents the limit of the applicability of the soft- 
pion dominance approximation as presented in this work. Description of processes involving 
larger momentum transfers would involve the introduction of the A excitation current, pion 
production processes as well as short ranged meson exchanges. There exists prescriptions 
to extend the present technique to incorporate these processes as shown in [ P^p3[ albeit in 



the non-relativistic SU (2) formalism. It would be interesting to extend these prescriptions 
to the present relativistic SU{3) formalism since it can trivially be used to estimate MEC 
corrections to any linear combination of SU{3) axial and vector currents shown in Eq. ( p.ll ). 
For example, the extended formalism may be used to investigate parity violating electron 
scattering processes which involves the interference of electromagnetic and weak neutral cur- 
rents. Other phenomena of interest to examine are neutral current scattering of solar and 
supernova neutrinos off deuterons [^| assuming finite strange quark form factors for the 
nucleons and inelastic nuclear transitions which can only take place in the presence finite 
strange quark form factors [^. The latter investigation would involve a relativistic finite 
nucleus neutrino scattering calculation which is the subject of the upcoming paper. 
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APPENDIX A: CONVENTIONS AND FORMULAE 

In this Appendix conventions and formulae used to compute semi-inclusive differential 
cross-sections for neutral and charged current neutrino-nucleus reactions in the RFC model 
are summarized for completeness. Because the technical details for including the Iplh 
and 2p2h MEC contributions to the hadronic tensor in the Fermi Gas models are given in 
detail in [^111) only the impulse approximation results are presented here. Although these 
references deal only with electron scattering, generalization to both charged and neutral 
current neutrino scattering is straightforward. Conventions of Halzen and Martin |58| are 
used throughout. 

As usual, it is assumed that the target nucleus consists of A nucleons with mass and 
that a single nucleon is ejected from the target due to the interaction with the probe. In 
what follows M is the mass of the free nucleon and the initial and final four momenta of the 
ejectile, = (£'/,Pj) and Pp = {Ef,Pf), are related to the four momentum transfer to 
the target k'^ by = — Pp. Incoming and outgoing lepton four momenta are denoted 
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— * — * 

by fcf = (ej, ki) and kj = (ej, kj), respectively. A positive energy nucleon with energy E in 
a plane wave state is given by 



V2EV 

where the 4-component spinor u{P, a), normalized to m^m = 2E, is 



(Al) 



u 



(P, a) = VE + M a-P 



(A2) 



.E + M> 



Here is the usual 2-component Pauli spinor with spin index a and the plane waves are 
normalized in a box of volume V such that 



V 



dr{rPw)Hpw = ^ 



(A3) 



The transition matrix element for the neutrino-nucleus scattering reaction is given in 
the usual way in the limit of very heavy vector boson masses [i.e. m|o iy± ^ 



Sfi = ig"^ Id'^x i<fy 



d^k 
(2^ 



(A4) 



where j^{x) and J^{y) are the leptonic and hadronic currents, respectively, and the coupling 
constant g is defined as g^ = Gp/^/2. The spatial parts of the leptonic and hadronic currents 
in the above transition matrix are 



fix 



Yj^ uikf, CT/)7^(1 =F l5)u{ku o^i) 



2V ie,ej) 
I (,i{Pi-PF)-y 



(A5) 
(A6) 



The minus and plus signs in Eq. ( |A5[) correspond to neutrino and anti-neutrino scattering, 
respectively. In the impulse approximation, the hadronic current operator J^{k) is given by 
either Eq. ( p.8|) for the neutral current or by Eq. ( p.9|) for the charged current, while the 
corresponding two-body MEC operators are discussed in Appendix B. The resulting cross 
section for free nucleons is given by 



(da 



Free 



6'^^\k^ - k^ + P^ - P^) a 



W^"" d^Ppd^kf (A7) 



where and cr^* assume the following form for neutral and charged current processes 

2 



a^' = 16e',cos'(^/2) 



47r 



a 



= l&k) 



An 



(A8) 
(A9) 
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with 9 being the lepton polar scattering angle. The quantity uo^yW^" appearing in Eq (|A7|) 
may be expressed as a linear combination of various nuclear structure functions Wj as 



(AlO) 



These structure functions are most conveniently expressed in a coordinate system defined 
by the following orthogonal unit vectors 



z = k 



n_L = 



kx Pp 
\kx Pf\ 

h± X k 
\h± X k I 



(All) 
(A12) 

(A13) 



In this coordinate system the spatial part of the hadronic current J, where = {p, J ), 
may be written as 

J = Jpii + J±h± + Jkk (A14) 
Then the structure functions for the semi-inclusive neutral current (z/, v'N) reaction are 



UJrpWT 

w± 



\p\' + 



UJ 



I T |2 



2Re {p*Jk 



uj±W± + uJi\W\\ 



\Ji 



Wtt = sm{2(f)F)Re (J^ jf^ 



LT 



Re 



(p - j^Jk) sin (j)pj\_ + cos (j)F J 



Wtt' = ImiJ\\J 



W. 



LT' 



Im 



(j) —Ji?j (^sin0j7j|| + cos </)irj| 

}\, 



(A15) 

(A16) 
(A17) 
(A18) 
(A19) 

(A20) 

(A21) 

(A22) 



where (l)p is the azimuthal angle of the ejected nucleon. The kinematical coefficients cuj 
corresponding to these structure functions follows 



= tan2(^/2) - 



1 + COs(20i;') 

1 — cos(20ir) 



.2\k\\ 

.2\k\\ 



UJtT 



(A23) 
(A24) 

(A25) 

(A26) 
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ujlt = -2 



\ 



^2 

tan(0/2) - ^ 



\k 



cott' = T2tan(^/2) 
ult' = ±2tan(e/2) 



\ 



tan(^/2) 



1.2 



\k 



(A27) 

(A28) 
(A29) 



In the expressions for the coefficients ujtt' and ujlt', the upper and lower signs correspond 
to neutrino and anti-neutrino scatterings, respectively. For the charged current reaction of 
type (uiJN) with lepton mass mi, each term in the product uj^^W^^, Eq. ( [A10[ ), may be 
written as 



4:eikf 
+ 



1 - kW 



m, 



\P\ 



2 I 2 

uj + rrii 



\k 



\Jk 



2uj\k 



Re {p*Jk 



(A30) 
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(A31) 
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-J^j J|| sin ^iT' + Jj_ cos ^i?) 
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(A32) 
(A33) 
(A34) 
(A35) 



In both the charged and neutral current processes the 4-momentum of the scattered 
lepton is not observed and it is thus necessary to integrate over the unobserved final lepton 
3-momentum kf in addition to the angular variables of the ejected nucleon. This leads to 
the following differential cross section for free nucleons 



\ dEp 



a^"/^* 27r 



Free 



Ae 



f 



\k 



(A36) 



In the relativistic FGM, non-interacting nucleons are ffiled up to the Fermi momentum kp 
and the initial momentum of the struck nucleon have to be averaged over the Fermi sphere. 
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In addition, the initial energy of this nucleon Ej is often reduced with respect to its free 
space value by the binding energy, expressed by an input parameter B 



Ej = {\Pi\'^ + M'^]^^^ - B 



= {\k\'^ + \PF\'^-2k-PF + M^y^^ -B (A37) 

The resulting differential cross sections per nucleon for neutral and charged current processes 
in the RFGM are 

V ^p^^ [4k].) Jo ^ ^ ' \k 



/ defd(cos6) — ujlWl + ujtWt + Wtt> (A38) 

Jo ' \k\ ^ i ^ ' 



Integrations over ej and cos^ in Eq. ( |A38|) are performed numerically. For further technical 
details see 



APPENDIX B: TWO BODY EXCHANGE CURRENT OPERATORS 

This Appendix presents the general form for the two-body soft-pion MEC operator, 
J^{k; Pi^i, Pj^2', Pp,!] Pf,2) EX, given in Eq. ( p.l6|) where the probing current, J^{k), may be 



any one of vector or axial SU{3) currents. Since the matrix element for the pion absorption 
is given by Eq. (|2.17|) , the quantity of interest here are the amplitudes for the soft-pion 



production shown in Eq. (|2.20|) . Therefore, only these amplitudes are shown explicitly here 
and the relevant full soft-pion MEC operators for electromagnetic, neutral and charged cur- 
rent scattering reactions may be constructed by taking the appropriate linear combinations 
of the operators given in this appendix, multiplied by the pion propagator and the nNN 
vertex function as shown in Eq. ( p.l6 ). 



Using the notation introduced in the text and in Appendix A, the amplitude for soft-pion 
production induced by the vector singlet current V° is found to be 



\un{7r^{q)N{PF)\V^{k)\N{Pj)) = 

{-i)u{PF, ctf) [VAl5lolf. + V^757m7o + V^757oS^ 

+V275S^7o + VIj.J:,] X^Pj, aj) (Bl) 

In the above expression a = 1,2,3 for pion production and the coefficients Vj multiplying 
the operators are defined as 

(B2) 
(B3) 
(B4) 
(B5) 
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9ttNN 

M 



(B6) 
(B7) 



with Ei = Y l-Pi P + M^. The vector singlet amphtude, Eqs. ( pTD is divergenceless when the 
nucleons are assumed to obey the free Dirac equation, leading to current conservation for 
vector singlet MEC. The corresponding vector octet soft-pion production amplitude is 

\im{7ryq)N{PF)\V'{k)\N{Pj)) = 



+ 4 



A", A' 



A^A^ 



V1757o7m + V^i757M7o + V^c757oS^ + ^^^75^^70 + ^^75^, 

(^a757o7m - '^1757^70 + 1^c757oS;. - 1/^75X^70 + 1/^757^ 

+Vo'j5k^)}u{Pi,cxi) (B8) 



where the index b runs from 1 to 8. The coefficients Vj are 
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I QirNN 
f Q-kNN 
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V El 

V M 

2. , G\{e 



A;2 




(B9) 
(BIO) 

(Bll) 
(B12) 

(B13) 

(B14) 

(B15) 



The first term in the coefficient Vq yields the "pion-in-flight" MEC with the pion form 
factor FTj{k'^Y. However, this "pion-in-flight" term is proportional to the four momentum 
transfer of the leptonic probe, /c^, and thus does not contribute to the neutral current total 
cross section. It also does not contribute to the electron scattering total cross section when 
the usual approximation of neglecting the electron mass is made. Terms proportional to 

also do not contribute to the (e, e') longitudinal and transverse response functions Rl 
and Rti when they are evaluated by using all the components of the electromagnetic four 
current as done in this paper. The coefficient Vq can be made to vanish by assuming a 
point like pion and using the pion pole dominance approximation from PCAC, Eq. ( p.ll ). 
For momentum transfers ranging from /5-decay to the quasi-elastic kinematics where the 
soft-pion MEC operators are applicable, this should be a good approximation. 

However, even when Vq is made to vanish the vector octet soft-pion production ampli- 
tude, Eq. ( pSD , is not divergenceless. The problem is the term with coefficient Vp which 
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has the structure of an axial current. The prescription used by Adler to make Eq. 
divergenceless and thus guarantee the conservation of vector currents involves adding ap 
propriate counter terms proportional to k^. In the present case, this counter term is ^^75^^ 
where Vq is given by 



(B16) 



This prescription is applicable to neutral current and electromagnetic reactions since terms 
proportioinal to fc^ do not contribute to the total cross sections. However, it leads to am- 
biguities for the charged current reactions where one of the leptons is massive and terms 
proportional to in the hadronic current must explicitly be taken into account in the deriva- 
tion of the reaction cross section. Nevertheless, since the generalized soft-pion dominance 
method of Chemtob and Rho ought to be applicable for all types of currents, it is desirable 
to insure the conservation of vector currents in this manner. In the present application, the 
contribution from the counter term VQ'j^k^ to the charged current cross section was found 
to be numerically very small (about 1% of the total MEC contribution). 

The matrix element of interest for the probing axial singlet current may be written as 



lim{n%q)N{P^)\A'^{k)\N{Pj)) = 

{+i)u{Pf, ap) Ua7o7m + ^b7/.7o + ^c^/. + ^2p7oM A"u(P7, ai) 
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(B18) 
(B19) 
(B20) 
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Similarly, the corresponding amplitude for the axial octet current is 
\iM'i^\q)N{P,)\Al{k)\N{P,)) = 

q— >U ^ 



{-^i)u{PF, CTp) 



with the following coefficients 



(A^7o7m - ^b7^7o + ^4^70^;. + A^Pi + ^2)^ + ^^7^) | m(P,, a,) (B22) 



A'a - +G\{e) (^) 



(B23) 
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(B28) 



(B27) 



(B24) 



(B25) 



(B26) 




(B29) 



Note that it is straightforward to generahzc these soft-pion production amphtudes to pseu- 
doscalar octet soft-meson production amphtudes where the SU{3) index a now takes on 
values from 1 to 8 and the ttNN couphng constant is replaced by the appropriate meson- 
nucleon coupling constants. However, in deriving the above amplitudes the Goldberger- 
Trieman relation was used to simplify the expressions and this algebraic manipulation is not 
necessarily justifiable for other pseudoscalar octet mesons under consideration. 
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FIGURE CAPTIONS 



FIG 1. a) Longitudinal, i?^, and b) transverse, response functions for the ^^C(e, e') 
reaction plotted against the energy transfer uj. The three momentum transfer involved in the 
inclusive reaction is |/c | = 400 MeV. A RFG model without any binding energy corrections 
[B = 0) and with a Fermi momentum of kp = 225 MeV was used to calculate the response 
functions . Experimental data, taken from [Q, are shown for reference but no attempt to 
fit the data by varying B and kp has been made. The dashed line represents the impulse 
approximation results while response functions obtained with soft-pion MEG corrections 
are shown with solid lines. Also shown in the figures using the dashed-dotted lines are the 
individual contributions from soft-pion two-body courrents. 



FIG 2 C{i'^- , jj )X total cross section obtained by folding the LSND neutrino flux 
3^ versus the Fermi momentum kp- The dashed curve is the impulse approximation result 



while the solid curve is obtained with the soft-pion MEG corrections using the RFG model 
of the nucleus with no binding energy correction {B = 0). For kp = 225 MeV, which is 
the usual value used for ^^C, the total cross-section is reduced from 24 to 22.7 (xlO~^°cm^) 
when corrected for two-body current effects. 



FIG. 3 a) ^'^C{u^- , ij,~)X total cross section as a function of incoming neutrino energy 
El, taking into account the Goulomb correction for the outgoing muon. As in Figure 2, the 
RFG model is used to model the nucleus with kp = 225 MeV and B = 0. The dashed 
line is the impulse approximation result while the solid line includes two-body soft-pion 
MEG effects. Note that in this particular case there is very little difference between the two 
results, b) Differential cross section da/dE^ for the same process obtained with and without 
the soft-pion MEG correction plotted against the outgoing muon kinetic energy E^. The 
results for this figure is obtained by folding the LSND neutrino energy distribution obtained 
from IBU . 



FIG 4 ^^C(z/, z/'p) differential cross section versus the kinetic energy of the ejected 
nucleon, Tp, for a) kp = 200 MeV, b) kp = 300 MeV and c) kp = 350 MeV. The incident 
neutrino energy is assumed to be 200 MeV while the values for the strange quark form factors 
used are F2 = —0.21 and G\ = —0.19. The long dashed curve is the impulse approximation 
result while the solid curves have been obtained with soft-pion MEG corrections. 



FIG. 5 ^^C(z/, z/'p) differential cross section versus the kinetic energy of the ejected 
nucleon, Tp, with {F2 = —0.21 and G\ = —0.19) and without finite strange quark form 
factors. The incident neutrino energy is assumed to be 200 MeV and the long dashed curve 
is the impulse approximation result while the solid curves have been obtained with soft-pion 
MEG corrections. 



FIG. 6 Various structure function contributions to the ^^C(z/, u'p) differential cross sec- 
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tion shown in Figure 5. Both the impulse approximation and MEC corrected results for the 
transverse, utWt, longitudinal, ujlWl, and transverse-transverse, utt'^tt', contributions 
are shown explicitly using the dashed and solid lines, respectively. Note that the two-body 
corrected contributions from the transverse and transverse-transverse contributions cancel 
leading to a small overall MEC correction to the impulse approximation. For anti-neutrino 
scattering the transverse-transverse contribution changes sign resulting in a smaller differ- 
ential cross section but in a larger MEC effect compared to neutrino scattering. 

FIG. 7 ^'^C{i>,i>'p) and ^'^C{u,ij'p) differential cross sections versus the kinetic energy 
of the ejected nucleon, Tp. The incident neutrino energy is 200 MeV and the nucleon is 
assumed to have strange quark form factors of F2 = —0.21 and G\ = —0.19. The long 
dashed curve is the impulse approximation result while the solid curves have been obtained 
with soft-pion MEC corrections. 

FIG. 8 a) Ratios of integrated proton-to-neutron quasi-elastic yield for the ^^C(z/, u'N) 
reaction as functions of for two values of strange magnetic form factor F2. In each case, 
the dashed line is the impulse approximation result while the solid line has been corrected 
for meson exchange currents. The incident neutrino energy is assumed to be 200 MeV for 
both cases and the range of integration was chosen to be 50 < Tp < 120 MeV to simulate 
the LSND experiment [^. b) Same as in a) but for anti-neutrino scattering. 
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